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Shape-shifting droplet networks 
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We consider a three-dimensional network of aqueous droplets joined by single lipid bilayers to 
form a cohesive, tissue-like material. The droplets in these networks can be programmed to have 
distinct osmolarities so that osmotic gradients generate internal stresses via local fluid flows to 
cause the network to change shape. We discover, using molecular dynamics simulations, a reversible 
folding-unfolding process by adding an osmotic interaction with the surrounding environment which 
necessarily evolves dynamically as the shape of the network changes. This discovery is the next 
important step towards osmotic robotics in this system. We also explore analytically and numerically 
how the networks become faceted via buckling and how quasi-one-dimensional networks become 
three-dimensional. 


In nature there are many biomaterials that are inter¬ 
nally programmed to morph into complex structures that 
actively and adaptively interact with the environment. 
Polypeptide chains fold into proteins and tubular lipid 
membranes self-assemble in branching networks to form 
the endoplasmic reticulum. Organisms as a whole de¬ 
velop shape by reorganizing the spatial distribution of 
their constituent cells in morphogenesis. 

These biological examples have inspired the develop¬ 
ment of programmable materials that controllably fold 
into designated structures. At the nanometer scale, the 
programmable chemistry of Watson-Crick base pairing 
allows DNA to self-assemble into a tetrahedron pQ. At 
the millimeter scale, efficient algorithms have been con¬ 
structed to generate self-folding three-dimensional poly- 
hedra from two-dimensional nets, driven by the mini¬ 
mization of the surface tension of liquid hinges that either 
rotate or fuse panels into place [2;. At the centimeter 
scale, researchers have created a self-folding robot that 
goes from flat to walking in several minutes without ex¬ 
ternal intervention [3]. 

A beautiful realization of these ideas comes from the 
Bayley group who print tens of thousands of micron-sized 
aqueous droplets each joined by single lipid bilayers lit 
8 to form a cohesive, tissue-like material (Fig. [l]). The 
droplets in these networks can be endowed with different 
osmolarities. The resultant osmotic pressure leads to lo¬ 
cal fluid flow from low to high concentration. This swells 
the high concentrations droplets and shrinks the low con¬ 
centration droplets, leading to internal stresses which dis¬ 
tort the shape of the network in specific ways depending 
on the initial geometry of the network and the concentra¬ 
tion differences. The initial droplet network experiments 
explored the formation of a hollow sphere from an ini¬ 
tial two-dimensional four-petal-shaped structure (see SM 
Movie-1), thus demonstrating the spontaneous assembly 
of three-dimensional shapes. 

To push the capability of this system into the realm 
of osmotic robotics, the droplet network must be able 
to make reversible shape changes. We investigate this 
possibility by turning to programmable interactions with 
the network’s environment. Past examples of a system’s 
environment driving shape change do exist. A strip of 
paper will, for instance, spontaneously curl up in your 


hand due to the interaction between the paper and the 
moist evaporative boundary layer of the hand [Q]. Bacil¬ 
lus spores respond to relative changes in environmental 
humidity with low humidity causing the spores to shrink 
and high humidity causing the spores to expand m- 
When these spores self-assemble into a dense monolayer 
sitting on a substrate, one can cycle the relative humidity 
of the environment to form an actuator. Osmotic actu¬ 
ation is used in plants, perhaps because it can generate 
a variety of plant movements, depending on the environ¬ 
mental conditions m, without consuming much power. 

Here we use predominantly computational methods 
to explore this plant analogue in the synthetic, micron¬ 
sized droplet network and show how shape shifting can 
be made reversible by exploiting the dynamically evolv¬ 
ing coupling to the environment, which here is simply 
the surrounding solvent. We illustrate reversibility with 
shape shifting from a four-petal configuration to a hol¬ 
low sphere transition and back. We also demonstrate 
the generalizability of this property by introducing first 
the formation of tetrahedral shells via folding and then a 
reversible tetrahedral shell folding-unfolding-folding pro¬ 
cess. Tetrahedra, the simplest and least symmetric class 
of the regular polyhedra, can serve as mesoscopic build¬ 
ing blocks for molecules and bulk materials with sp 3 -like 
directional bonding [12, 13. Our reversibility findings 
open the door to osmotically-driven small scale robotics. 

Finally, to better determine the range of shape forma¬ 
tion in these droplet networks, we revisit the formation 
of rings studied in Ref. |4 and identify a new buckling 
transition from a round to a polyhedral ring after ring 
closure. This configuration also has the advantage of 
being analytically tractable. We also explore spiral for¬ 
mation as they provide an interesting example of an ini¬ 
tial quasi-one-dimensional structure generating a three- 
dimensional structure to expand the self-assembly capa¬ 
bility of the system. 

Model- We model the droplet network using molecu¬ 
lar dynamics in three-dimensional Euclidean space. Each 
droplet i is treated as a point mass m^, with an associ¬ 
ated radius Ri and osmolarity (defined as the number of 
osmoles (Osm) of solute per litre (L) of solution) Q. As 
described in Ref. [4], each droplet interacts with neigh¬ 
boring droplets via an elastic interaction and an osmotic 
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FIG. 1: (Color online) Schematic image showing droplets of 
different aqueous solutions printed into a solution of lipids in 
oil. The droplets acquire a lipid monolayer and form bilayers 
with droplets in the developing network. 


interaction. The elastic interaction potential of a pair of 
droplets i and j via the bilayer attaching the two droplets 
is written as 

„ f \k [rij - l (Ri + Rj )] 2 for rij < Ri + Rj 

Eij — \ ^ 

[ 0 for > Ri + Rj 

(i) 

where k is the spring constant, is the distance between 
a pair of droplets i and j, and l = 0.8 represents the 
change of the equilibrium length due to the deformation 
of two droplets when they are fused [4]. A damping force 
for each droplet proportional to its velocity is included 
as demanded by Stokesian flow. The exchange of water 
between droplets of different osmolarities is described by 
Fick’s first law as 


Jij = AijD (Cj — Ci ), ( 2 ) 


where Jij is the volume of water transferred per unit 
time from a droplet i with osmolarity Ci to a connected 
droplet j with osmolarity Cj , Aij is the common inter fa¬ 
cial area of two connected droplets i and j and D is an ef¬ 
fective permeability coefficient taken to be constant. The 
flow of water between two connected droplets changes 
the size of each droplet. To simulate the dynamics of 
the droplet network, the net force on each droplet Fi is 
described by 
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where Vi is the velocity of droplet i and 7 is the damping 
coefficient. In experiments, the radius of the droplet is of 
order tens of microns, the spring constant is of order a few 
milliNewtons per meter and the time scales for mechani¬ 
cal relaxation and water transfer are seconds and tens of 
minutes, respectively. We invoke this approach because 
the time scale for mechanical equilibrium is much faster 


than the time scale for water transfer. More simulation 
details may be found in SM II. 

Reversible Folding- Let us start with the four-petal¬ 
shaped structure studied experimentally in Ref. 0- By 
adding an osmotic interaction with the environment 
we can realize a reversible folding-unfolding process, as 
shown in Fig. [ 2 ] (see SM Movie-5). To generate reversible 
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FIG. 2: (Color online) (a)-(c) Three snapshots from the re¬ 
versible folding process. See SM Movie-5 for the full sequence 
of snapshots. 

folding we place part of the folded droplet network into a 
medium with higher osmolarity C m (C m > C 72 ) so that 
water will flow from the droplets on the bottom layer of 
the “flower” to the medium. More precisely, the bottom 
layer droplets located below a horizontal x — y plane at 
z = h are exposed to the higher osmolarity medium. As 
the flower folds, an increasing number of bottom layer 
droplets naturally become exposed to the medium. They 
therefore lose water, to the surrounding medium, and 
start to unfold again. The overall volume of the unfolded 
flower also drops, preventing a complete reversal. As the 
bottom droplets continue losing water to the medium, the 
top droplets are also losing water to the bottom droplets. 

The detailed evolution depends on the the osmolar¬ 
ity difference and the total contact area through which 
the water is passing. Lower ft, which means the smaller 
the contact area between the bottom droplets and the 
medium, and the smaller C m , which means a smaller os¬ 
molarity difference, will lead to a slower rate of water 
transfer between the bottom droplets and the medium. 
For top layer osmolarity C top = 0.1 and bottom layer os¬ 
molarity Cbottom = 1 .0, one finds that the top layer loses 
more water to the bottom layer than the bottom layer 
itself loses to the surrounding medium. This eventually 
reverses the unfolding and the structure starts to fold 
once again. As discussed in SM II, our simulation uses 
the same algorithm as before, other than a modification 
of Fick’s first law (Eq. ( 2 )) resulting from the separation 
of time scales between mechanical relaxation and water 
transfer. 

To characterize the reversibility of the system, we de¬ 
fine reversibility as 

P = (E m ax-before T Dmax-after c 2D rn i n ) /iF norm , (4) 

where D m i n is the minimum depth of unfolding, 
E max-before is the maximum depth of the flower before 
the “flower” unfolds, D max _after is the maximum depth 
of the flower after it reaches a minimum and refolds, and 
Enorm is a normalization constant. For p = 1, the sys¬ 
tem is fully reversible in the sense that the flower opens 
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out fully before folding back up again, while for p = 0 
there is no reversal. 

Fig-! shows the dependence of p on the horizontal 
plane’s z = h coordinate and the media’s osmolarity C m . 
As the horizontal plane rises, the reversibility increases 
since there are more bottom droplets in contact with the 
medium. The reversibility also increases with a rise in 
the media’s osmolarity C m . 


FIG. 3: (Color online) Reversibility dependence on the z co¬ 
ordinate of the horizontal plane and the osmolarity of the 
surrounding medium. 

We can generalize this notion of reversibility to tetra¬ 
hedral shells folding and unfolding. Before addressing the 
reversibility, we first demonstrate how a tetrahedral shell 
can form. We choose as initial droplet network a central 
triangle connected to three other triangles by hinges on 
each of its sides. This is an example of a two-dimensional 
net of polygons and hinges (creases in origami) that can 
be folded into three-dimensional structures. We use four 
layers of droplets, all of osmolarity Ci, to create suffi¬ 
ciently rigid faces. For each hinge, the top two layers of 
droplets have osmolarity C \, while the bottom two layers 
droplets have osmolarity C^> C\. The flow of water to 
the outer layers causes the hinges to bend upward, clos¬ 
ing all faces into a tetrahedral shell (see Fig. [4] and SM 
Movie-4). The osmolarity difference must be tuned to 
achieve tetrahedral folding. If C\ = 1.0 and C 2 < 4.0, 
the osmolarity difference is too small to drive tetrahedral 
closure. This threshold can be approximated analytically. 
See SM III for details. The critical osmolarity difference 
may be lowered or raised by using wider or thinner hinges 
respectively. 

Since this structure has yet to be explored experimen¬ 
tally, we also test for robustness of the tetrahedral shell 
formation to disorder. We do so by randomly diluting 
droplets in the initial state. For 20 percent dilution and 
lower, tetrahedral shell formation persists. This percent¬ 
age can be increased as the width of the hinges increases. 
Moreover, should we allow for one dilution percentage in 
the hinges and a second, independent dilution percentage 
in the triangles, we expect the percentage of dilution in 
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FIG. 4: (Color online) (a)-(d) Forming a tetrahedral shell 
(Ci = 1.0 and C 2 = 5.0). The color represents osmolarity 
with blue/red representing low/high osmolarity with white 
intermediate. 


the triangles at which the tetrahedral shell formation is 
no longer robust to be higher. 

And now for the reversibility capability of the tetrahe¬ 
dral shell. See SM Movie-6. If we assume that only the 
bottom layer droplets of the hinge transfers water with 
the solvent, there is a smaller window of the osmolarity 
of the solvent C m than the four-petal-shaped structure 
for the reversible folding to occur. In particular, the 
relation between reversibility capability and C m is 
nonmonotonic, e.g. for small and large values of the 
osmolarity of solvent, the folding is not reversible. For 
the larger values of C m , the triangles in the tetra¬ 
hedral shell will not be flush as the volume of the 
hinge decreases faster than the triangle. Moreover, the 
reversibility depends on the height of the surface of the 
solvent strongly than the four-petal-shaped structure be¬ 
cause the hinge is shorter than each petal of the “flower”. 

Rings and buckling- The buckling of a circular elastic 
ring subject to an external radial pressure has been ex¬ 
tensively studied in applied mechanics [14]. This buckling 
can presumably occur due to internal pressures as well. 
To study this possibility in the context of droplet net¬ 
works we evolve from an initial configuration consisting 
of two rows with different osmolarities (see SM Movie-2). 
FIG. [5] shows ring closure for N = 38 total droplets. The 
initial osmolarity of the top row is C\ =0.1 (blue) and 
that of the bottom row is C 2 = 1.0 (red). After each row 
closes to form a ring there is still an osmolarity difference 
between the outer and inner rings, as can be seen from 
Fig{5|c). This residual osmolarity mismatch is followed 
by a ring buckling transition, as shown in Figj5je). 

The final shape depends on the osmolarity difference 
and the number of droplets. In Fig. [6| we sketch the 
phase diagram for ring closure and buckling as a function 
of the top row osmolarity C\ and the total number of 
droplets N. The bottom row osmolarity is fixed at C 2 = 
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FIG. 5: (Color online) (a)-(e) Buckling of a ring with N = 38 
total droplets. The initial osmolarities of the two rows are 
Ci =0.1 (blue) and C 2 = 1.0 (red), respectively. Snapshots 
were generated using the Visual Molecular Dynamics (VMD) 
package m and rendered using the Tachyon ray tracer m ■ 



FIG. 6: (Color online) Plot of the ring buckling phase diagram 
as a function of top row osmolarity C\ and total number of 
droplets N. The bottom row osmolarity is fixed at C 2 = 1.0. 
The symbols, obtained from simulation, should be compared 
with lines showing analytical results. 


1.0. There are three phases. For a given value of TV, 
rings do not form at all until the osmolarity difference 
AC = C 2 — Ci exceeds a threshold. The threshold value 
is computed analytically in SM IV. For a range of AC 
one then finds smooth rings. For yet larger AC the closed 
ring buckles. An upper bound for the buckling transition 
can also be obtained analytically (see SM IV ). The larger 
AT, the easier it is to form a ring and so the smaller is 
the threshold osmolarity difference for ring closure and 
subsequent buckling. 


FIG. 7: (Color online) Typical initial (left) and final (right) 
configurations for the formation of spirals. 


Spiral formation- Finally, to explore the possibility of 
generating three-dimensional structures from quasi-one- 
dimensional structures, we explore spiral formation. We 
choose as initial state a 200 x 5 x 5 rectangular slab (see 
Fig.[7]and SM Movie-3). Each droplet can be indexed by 
integer orthogonal coordinates /, m, n, with 0 < m < 200, 
0 < n < 5 and 0 < l < 5. The initial osmolarity of each 
droplet is set to be 

Cmni = (1.0 - K n n ) ( 1.0 + Ki l (m ~^° 0) ) , (5) 

where K n and Ki are free parameters. We explore how 
spiral formation process depends on FQ and K n for fixed 
slab size. 



FIG. 8: (Color online) (a) The radius of the spiral as a func¬ 
tion of K n for different Kis. (b) The pitch of the spiral as a 
function of K n for different Kis. 

Fig. J8^a) shows how the radius of curvature of the spi¬ 
ral depends on FQ and K n . K n determines how efficiently 
the slab folds up in the m — n plane. As K n increases the 
rod folds up more effectively and the radius of curvature 
decreases. A positive value of FQ will make each part of 
the slab fold up differently in the l direction. If FQ =0, 
for example, the two ends of the slab will fold and then 
meet each other, leading to a ring. Note that the radius 
of the spiral also decreases with increasing FQ, making 
folding more efficient. 

Fig. |8jb) plots the pitch of the spiral as a function of 
Ki and K n . The pitch depends mainly on FF n , with more 
efficient folding occurring at large K n . There is a slight 
dependence on FQ. 

Discussion — We have demonstrated for the first time 
how to, in principle, reverse the large-scale shape change 
in these droplet networks with the four-petaled flower 
folding, unfolding, and folding back again. We have 
also demonstrated the generalizability of this reversibility 
property to the folding-unfolding-refolding of tetrahedral 
shells. Forming a general shape from a net with folds 
in general (and the subsequent unfolding and re-folding) 
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requires simultaneous fold movements. This should be 
possible to achieve with the right set of osmolarity param¬ 
eters. The analytical analysis of the hinge in SM III also 
paves the way for analytically predictive estimates of the 
formation of more general shapes. Our folding-unfolding¬ 
refolding processes capitalizes on the interaction between 
part of the droplet network and the surrounding medium 
(environment) and relies on the medium having an osmo¬ 
larity larger than the osmolarities in the droplet network. 
Reversibility is a first step towards osmotic robotics, 
namely a robotic gripper. The reversibility of our model 
is limited, however, because (1) the droplet network folds, 
unfolds, and folds again only once and (2) the volume of 
the final state of the folded shell (spherical or tetrahedral) 
is smaller than if it had just simply folded. 

To address the limited reversibility of our model, a 
recent study extends osmosis to active solutes con¬ 
taining, for example, self-propelled colloids [151121] or hot 
nanoparticles [22j i2j3j. The study finds that active solute 
activity increases the osmotic pressure and can also expel 
solvent from the solution. By using active solutes, the sol¬ 
vent can be controlled to flow from the higher osmolarity 
to lower, which is reversed from the usual passive solute 


situation. The use of reverse osmosis via active solutes 
could make the folding-unfolding process completely re¬ 
versible. Further study, however, is needed to assess the 
feasibility of active solutes to achieve such a goal. 

We have also predicted the buckling of rings and the 
formation of spirals in these droplet networks driven by 
osmolarity gradients. For ring formation, we find a sub¬ 
sequent buckling transition for large enough osmolarity 
difference. For spiral formation, we determine how the 
size of the gradients in each of the three dimensions af¬ 
fects both the radius and pitch of the spiral. Both the 
programmable buckling and the quasi-one-dimensional 
to three-dimensional structure formation enhance the 
programmable shape-shifting capability of the network. 
Moreover, robots could use buckling to generate sudden 
changes of geometry and spiral formation to enhance ma¬ 
neuverability. 

The authors acknowledge useful comments by M. C. 
Marchetti on an earlier draft of the manuscript. MB 
thanks L. Mahadevan for stimulating discussions. This 
research was supported by the Soft Matter Program at 
Syracuse University. 
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Supplementary Material for “Shape-shifting droplet networks” 


I. SELECTED MOVIES 

1) A four pedaled structure folds into a flower shaped hollow sphere: 

Movie- l-folding_flower. mp4 

2) Ring formation and subsequent buckling: 

Movie-2-folding_ring.mp4 

3) A 200x5x5 rectangular slab folds into a Spiral: 

Movie-3-folding_spiral.mp4 

4) Tetrahetron formation: 

Movie-4- folding_tetr ahetron. mp4 

5) By adding an osmotic interaction with the environment, a folded flower shaped structure can 
unfold and then fold again: 

Movie-5-folding_unfolding_refolding.mp4 

6) A tetrahedron can also fold, unfold, and fold again: 

Movie- 6 -folding_unfolding_refolding_tetrahedron.mp4 


II. FURTHER SIMULATION DETAILS 


To simulate the dynamics of the droplet network, the net force on each droplet Fi is described by 


3 = -£ 


<ij> 


dEij 


• 7 Vi = m 


cPfi 

dt 2 


(SI) 


where Vi is the velocity of droplet i and 7 is the damping coefficient. Following Ref. [4], we assume m, &, and D are 
the same for all droplets. Because the time scale for mechanical relaxation is much faster than the time for water 
transfer (seconds compared to tens of minutes), any global shape change is in mechanical quasi-equilibrium. Thus 
the simplifying assumption of identical and constant m, fc, and D for all droplets should not affect our results in any 
significant way. We chose, in simulation units, D = 0.002, k = 10 3 , and m = 0.2. The value of the damping coefficient, 
7 = 1 . 1 , was chosen so that there are no oscillations between the droplets when they bind as observed experimentally. 
The osmolarities (C*) were chosen to give gradients similar to those for the flower-closing experiment. Length, time, 
and mass scales in simulation units can be converted to microns, seconds, and grams by matching to experiment. 

Implementing different initial osmolarity gradients and droplet configurations yields folding into a variety of impor¬ 
tant structures such as rings, spirals, and tetrahedra. To search for such structures, the droplets are initially positioned 
in hexagonal closed-packed arrangements with a common initial volume Vo = 10 (Vo = IttRq, with Ro the initial 
radius) and equilibrated first without water exchange, after which osmosis is switched on. Any two droplets in con¬ 
tact are then connected via the elastic interaction. Once the elastic interaction is established, water is exchanged via 
Fick’s first law and the radius of each droplet is updated accordingly (we assume the droplets always keep a spherical 
shape). The inter facial area Aij in Eq. | 2 | in the main text is approximated as A{j = ttR 2 ^ , where Rij = min {Ri, Rj}- 
The position of the center of each droplet is then updated using a fourth-order Runge-Kutta scheme to obtain the 
position of the centers at the subsequent time step with At = 0.01 in simulation units. We have checked that our 
simulation results are robust to making the time step as small as At = 0.001 and as large as At = 0 . 02 . As for the 
computational cost involved, a 4372 particle simulation for tetrahedron formation over 250,0000 simulation time steps 
took approximately 3 hours and 19 minutes on a computer with 2 quad core 2 . 66 GHz processors. 

To model the droplet network osmotically interacting with the solvent in a controlled way, we place the folded 
droplet network in a medium with osmolarity C m exceeding any individual osmolarity of the droplets. The surrounding 
medium can only interact with part of the droplet network, as described above. We did not simulate the medium 
explicitly and neglected any mechanical response of the medium. More simulation details and the experimental 
values of parameters can be found in Ref. [4]. 
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III. DETAILS OF AN ANALYTICAL ANALYSIS OF EACH HINGE’S FOLDING IN THE FORMATION 

OF A TETRAHEDRON 


While each hinge in the tetrahedron formation is a 3D structure of 4 layers hexagonal closest-packed lattice, here 
we analyze a 2D structure with only two rows of droplets as an approximation. Each row contains N droplets. The 
initial radius of each droplet is Ri and the final equilibrated radius of each droplet is Rif for the top row and R^f 
for the bottom row. The initial osmolarity is Cu for the top row, C<ii for the bottom row, and the final equilibrated 
osmolarity is Cf for both rows. The amount of solute in each droplet is conserved. This implies 


4 o 4 o 

Cu-nRi 3 = Cf-nRif 3 , 

4 , 4 , 

C^nRi 3 = C f -nR 2f 3 . 

The total amount of solution is also conserved: 


This yields 



fKR lf 3 + ^7Ti?2/ 3 . 


Rif 


i?2/ 


2 Cu 


Cl i + C^i 


( 2 C 2i \ 

\Cli + C^i) 



(S2) 


(S3) 


(S4) 


To fold this hinge by an angle #, note that each row of droplets must have the same center 


of curvature, labelled 



as point O in Fig. SI. The distance between two neighboring droplets in the top (inner) row is 2/ithj, the distance 
between two neighboring droplets, one from the top row and the other one from the bottom row, is l (Rif + i? 2 /), 
and the distance between two neighboring droplets from the bottom row is 2 /R 2/5 where / = 0.8 as shown in Fig. SI. 
In A OAB, we have 


\OA\ = ( IRif) / sum 
\OB\ = (/i? 2 /) / sum 
\AB\=l(R lf + R 2f ) 

\AB\ 2 = \OA\ 2 + \OB\ 2 - 2\OA\\OB \cos a 
This yields for onset condition of folding this hinge by an angle 6 = ( 2N — 1) a 

+ R\ s - 2R lf R 2f cos (j^zr) 


(Rif + R 2 f) — 


R\ f 


sin 


\2N-1 J 


(S5) 


(S6) 
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By combining Eqs. |S4| and S6 we can numerically solve for C 2 { if given the value of N. For example, when N = 4 
and Cu = 1.0, we plot the onset value of C 2 i to fold the hinge by an angle 0 as shown in Fig. S2. To form the 



FIG. S2: The folding angle of a hinge 6 as a function of the onset value of C 2 i, when Cu = 1.0 and N = 


4. 


tetrahedron, it requires 0 = 109.47°, then the onset value of C 2 { is 4.20. 


IV. DETAILS OF ANALYTICAL ANALYSIS OF RING FORMATION AND SUBSEQUENT BUCKLING 


We consider the ring formation of two rows with N total droplets, as shown in Fig. S3. In this situation, similar 

[S4l 


to hinge folding, Eqs. 


and 


condition of ring formation is 


S5 


are valid for the equilibrated configuration, with a = WJ2 in E( l- 


S5l 


And the onset 


(Rif + Rif) — 


Rh 


R\f 2RifR 2 f cos ^ (n/ 2 ) ) 


sin 


( to ) 


(S7) 


By combining Eqs. [S4|and[S7| if given the values of C 2 i and TV, we can solve for Cu numerically. 

To determine the parameters for the onset of ring buckling where the inner droplets become separated, note that 
the distance between two neighboring droplets from the top(inner) row is 2 Rif, the distance between two neighboring 
droplets one from the top row and the other one from the bottom row is l\ 2 ( Rif + R 2 f), and the distance between 
two neighboring droplets from the bottom row is l 22 ( 2R 2 f ), as shown in Fig. S4. For a droplet from the inner row, 
the force exerted by each of the two neighboring droplets from the same row is \Fu\ = (1.0 — l) (2i7ij), and the force 
exerted by each of the two neighboring droplets from the outer row is \F 2 i\ = (/12 — 0 {Ri / + R 2 /), where l = 0.8. 
Those four forces balance each other, yielding 


2|p ”' si "((W2))= 2 Al 


\J\l\2 (Rif + R 2 f)] 2 - (^ 22 # 2/) 2 

I 12 (Rif + R 2 f) 


(S8) 


so that 


/ 7r \ a/[/i 2 (Rif + R 2 f )] 2 — {I22R2f ) 2 

2 (1.0 - 0 (2Hi/)sm (In - l) (R„ + R„) * - lMf + Raf) - 


(S9) 


For a droplet from the outer row, the force exerted by each of the two neighboring droplets from the same row 
is \F 22 \ = (l — l 22 ) (2R 2 f) and the force exerted by each of the two neighboring droplets from the inner row is 

















FIG. S3: Schematic figure for the onset of ring formation. 


\Fl2\ 


(l \2 — l ) (Rif + where l = 0.8. Those four forces balance each other, yielding 


so that 


2 |F 22 | sin 



2|Fi 2 | 


\/[h2 (J?i/ + J? 2/ )] 2 -(J?i/) 2 

I12 (Rif + R-2 /) 


/ 7T \ \f [^12 (Rif + -^2/)] 2 — (Rif) 

2 (I - W (m .,)*. ( mj ) - 2 fe - 0 (a, + *,) v Ty^TTiwI - 

We also have 

l 0A l = ( E ^ /sin (^ 2 j) 

|OB| = (/ 22J R 2/ )/sin(^) 

= Z 1 2 (# 1 / + # 2 /) 

\AB\ 2 = \OA\ 2 + |OB| 2 - 2|OA||OB| cos (^y) 

The condition for the onset of ring buckling is thus 

„ o -^1/ + ^22-^2/ — 2i 22 i?i/i? 2 y COS f (jv/ 2 ) 1 

4 (iZ,, + i? 2/ ) 2 = —-- 7 -y- ^ 121 . 


(S 10 ) 


(Sll) 


(S 12 ) 


(S13) 


By combining Eqs. |S4[ S9 Sll and S13 we can numerically solve for Z 12 , Z 22 and Cu for given values of N and C 2( . 
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FIG. S4: Schematic figure for the onset of ring buckling. 





